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Abstract 

We study an analytically solvable model for decoherence of a two spin system embedded in a 
large spin environment. As a measure of entanglement, we evaluate the concurrence for the Bell 
states (Einstein-Podolsky- Rosen pairs). We find that while for two separate spin baths all four 
Bell states lose their coherence with the same time dependence, for a common spin bath, two of 
the states decay faster than the others. We explain this difference by the relative orientation of the 
individual spins in the pair. We also examine how the Bell inequality is violated in the coherent 
regime. Both for one bath and two bath cases, we find that while two of the Bell states always 
obey the inequality, the other two initially violates the inequality at early times. 
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Entanglement, nonlocal quantum correlations between subsystems, is not only one of 
the basic concepts in quantum mechanics y| but also central to quantum computation and 
quantum information |2|. Decoherence, loss of phase relations between the states, is essential 
in understanding how a quantum system becomes effectively classical [3]. Therefore, how 
an entangled system undergoes decoherence or how the entanglement changes as a result of 
interaction with the environment is an important issue and for two entangled spins subject 
to quantum noise created by a bosonic bath the problem has already been studied jj, 5[. 

In this work, we concentrate on decoherence of two spins as a result of an interaction 
with a spin bath. This problem is closely related to electron spin dynamics, due to hyperfme 
interaction with surrounding nuclear spins, in quantum dots |6(. Decoherence of various sys- 
tems, including superconducting quantum interference devices (SQUIDs) coupled to nuclear 
and paramagnetic spins, can be described by similar models [7]. Many spin systems can ex- 
hibit interesting behaviors including parity dependent decoherence where some nondiagonal 
elements of the density matrix survive the initial decay of other entries due to environmental 
spins QQ. For the centra, spin mode., wh.eh describes agnized spin coup.ed to a spin 
bath, the quasiclassical equations of motion are integrable JlJjJ. 

Quantification of entanglement is a major challenge in quantum information theory. A 
well known measure for a pure state of a pair of quantum systems is the von Neumann 
entropy or equivalently the Shannon entropy of the squares of the Schmidt coefficients jll|. 
The entropy of the partial density matrix, which is obtained by tracing out one of the 
members from the total density matrix, can be used to parametrize the entanglement. For a 
pair of binary quantum objects (qubits) an alternative parameter is the concurrence which is 
related to the von Neumann entropy bijectively Q|. To quantify the entanglement between 
the two spins, we are going to use the concurrence because of its mathematical simplicity. 
Our main results related to entanglement will turn out to be independent of the choice of 
the measure. 

Decoherence of the two spins can be viewed as a generation of entanglement between the 
pair and the spin bath (or baths) and hence any measure of the entanglement can also be 
used to parametrize the decoherence. As the members of the pair lose their entanglement 
with each other, they start to entangle with the bath spins. What we are going to evaluate 
is the concurrence corresponding to the entanglement between the two partners. 

Our aim is to understand how two entangled spins lose their correlation due to other 
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spins interacting with them. For this purpose we start with a very simple model where we 
can observe decoherence effects. The model Hamiltonian 

H = c u hu lk ai kz + c 2z Y hu 2 k°2kz (1) 

k=l k=l 

describes two central spins, with z— component operators C\ z and C2 Z , coupled to bath spins 
represented by <r n kz, where n = 1, 2 labels the baths and k = 1,2, 3, N n labels the individ- 
ual spins. All spins are assumed to be 1/2 and c\ z , ci z , and a nkz denote the corresponding 
Pauli matrices. Hamiltonian (|T|) is a simple two spin generalization of the model proposed 
by Zurek to study decoherence in spin systems ^|. First, we are going to consider two dif- 
ferent spin baths where each spin couples only one of them. Later, we are going to examine 
how our results change when the pair interacts with a single bath. Since the Hamiltonian 
involves only the z— components of the spins, it can also be used to study decoherence 
of other two-state systems. 

We are going to assume that at t — 0, the central spins are not entangled to the spin 
baths so that the state is in the product form |^(0)) = |^ c (0))|^ CT i(0))|^ o - 2 (0)) where 

l*c(0)) = (on | TT) + a u \ tl) + o iT | IT) + o u | ||» (2) 

with obvious notation for the two spins and 

l*«m(0)) = (g) Kfc| Tnfc) + Pnk\ ink)) (3) 
k=l 

where | ] n k) and | [ n k) are eigenstates of a nkz with eigenvalues +1 and -1, respectively, 
and \a n k\ 2 + \/3 n k\ 2 — 1- At later times, the state is no more in the product form due to 
entanglement of the pair with environmental spins but instead it is given by 

|*(t)) = (a tT I TT)|*al(+t)>|* ff2 (+t)> 

+a u | TJ.>l*«rt(+*)>l*aa(-*)> 

+a u | It)|*«xi(-0)l*<x2(+t)) 

+au I U)\**i(-t))\9o2{-t))) 



(4) 



where 

\^*n(t)) = (g)(a nfc e-^ *| Uk) + Pnke^l U»- (5) 

We are going to see that it is the randomness of the interaction strengths and the expansion 
coefficients that will lead to decoherence of the pair. 



The total, central spin and the baths, density matrix is given by pit) = \^(t)){9(t)\ but 
what we are interested in is the reduced density matrix which is obtained from the former 
by tracing out the bath degrees of freedom as p c (t) = Tr a p{t) where subscript a means that 
trace is evaluated by summing over all possible nk states. We can write the resulting density 
matrix in the product basis {| If), | fl), | IT), | II)} as 



' |a TT | 2 an a U r2 a TT a |T ri a TT a || r i r 2 
a| T a n r* |a n | 2 a^J^r* a n aj j r 1 
a| T a iT r* a^a^r*^ |a iT | 2 %r a u r 2 

M 2 



\ 



(6) 



V«TT a H r i r 2 a n a ll r i a h a ll r 2 
where * means complex conjugation and decoherence factors ri(t) and ^(t) are given by 



r n (t) 



n 

k=l 



; nk i 2 e i2uJnkt 



(7) 



In general both expansion coefficients (3 n k and interaction strengths uo n k are random. 
If the bath spins point randomly at t = we can write the expansion coefficients as 
ctnk = cos(6'„,fc/2)e _ ^" fc//2 and (3 n k = sin(6' n fc/2)e^ nfc/ ' 2 , where n k and 4> n k are spherical polar 
coordinates determining the direction of the spins, and we assume that the angles 9 n k and 
(p n k have uniform distributions in the intervals [0, ir] and [0,27r], respectively. It is possible 
to show that for sufficiently large N n values |r n (£)| exhibits a Gaussian time dependence 



rather than exponential 



14 1 . In our case 



16 ^2 I a nfc| 2 |/?nfc| 2 |^nfc| 2 - 



(8) 



We are going to obtain several coherence factors given by expressions similar to Eq. ((7j). We 
first note that the larger the interaction strengths |u; n &|, the faster the decay. Secondly, for 
a given set of {co n k}, the fastest decoherence is attained when \a n k\ and \(3 n k\ become equal. 

To evaluate the concurrence [l^j], we need to find the time- reversed or spin-flipped density 
matrix p c which is given by 



Pc 



0-,, 



Vy)p* c ((Ty <g> CT y ). 



(9) 



Here a y is the Pauli spin matrix and <g) stands for the Kronecker product, and p* c is obtained 
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from p c via complex conjugation. We can write the spin- flipped density matrix as 



/ 



Pc = 



ku\ 2 



-a iT a u r 2 -a n a uri a^a^r^ 

htl 2 



\ 



- a h a U r i a h a U r *i r 2 \a n \ 2 -a n a* n r 2 

|2 



(10) 



\atT°H r i r 2 -«TT a iT r i ~ a u a U r 2 l a ttl 

The final step in evaluation of the concurrence C is to find the four eigenvalues {Aj} of the 
product matrix pp c in the decreasing order so that 

(11) 



C = max{0, y Ai — y A 2 — y A 3 — y A 4 }. 

We are going to evaluate the above expression for the Bell states (Einstein-Rosen-Podolsky 
pairs) 

u \ _ im+iu) 



ITD+IIT) 

p \ _ im-iu) 

. [UHiD 



e 2 



e 4 



(12) 



As we are going to see, the Bell states have the property that the concurrence is the same 
for all of them. In fact, any other basis obtained from the Bell states by replacing the 
coefficients ±l/\/2 with e ie /y/2 (6 being a real number) has the same property. 

For two baths, the concurrence, which is the same for all of the Bell states, turns out to 

be 

C=\n\\r 2 \. (13) 

Since, r^O) = r 2 (0) = 1, the concurrence is also unity at t — 0. On the other hand, both r 1 
and r 2 , and hence the concurrence, decay with time and vanish. For the special case where 
only one of the spins, say the first one, interacts with a spin bath so that r 2 (t) = r 2 (0) = 1, 
we still observe a decay in the concurrence. This is an expected result because the entangled 
pair must be treated as a single system rather than individual spins. 

We next consider the case where both spins undergo decoherence due to interaction with 
the same spin bath so that the Hamiltonian becomes 



N 



H = Tl^itOikCiz + U} 2 kC2z)Vkz 
k=l 



(14) 
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This time the state at t = is given by |^(0)) = |* c (0)) ^(O)) where 



N 

1^(0)) = ®(a*| U)+Pk\ U)). 
fc=l 

Similar to the two bath case, we can write the density matrix as 

1 ' ' 12 an a u r 2 a TT a |T ri a TT a || r i2 ^ 

2 



(15) 



Pc 



a n a u r 12 a n a u n 

2 



|a n | 2 
«TT a Tl r 2 \ a U 

\a*K a n r i2 a n a ii r * a lt a ii r 2 l a u 

where decoherence factors are given by 

N 

r n (t) = II (K| 2 e-* 2 ^ + \p k \ 2 e^ 



(16) 



k=l 



and 



N 



12 



(i)=Il Kfe 



|2 -i2(w lfc ±w 2fe )f 



2„i2(cj lfe ±<j 2 fe)t 



(17) 



(18) 



k=l 



As we have discussed in the paragraph after Eq. (jSJ), the larger the interaction strengths, 
the faster the decay. Therefore, we should compare uo\k + 002k with uo\k — uJ2k- If all of the 
interaction constants u> n k have the same sign, r ] l 2 (t) goes to zero faster than rj^(t). For the 
special case, uik = u 2 k for all k, rf 2 (t) does not decay at all but remains constant. 

After finding the spin-flipped density matrix p c and eigenvalues of the product pp c , we 
can evaluate the concurrence for each of the Bell states. In single bath case the Bell states 
exhibit different decay rates with the concurrence expressions 



C\ — C3 

C2 = C4 



121 



121 



(19) 



Although we have obtained this two by two grouping of the Bell states in terms of the 
concurrence, any other measure, like the von Neumann entropy, which depends upon the 
eigenvalues of the density matrices will yield the same result. Equations ()13|) and (|T9~|) 
show that the concurrence, which is a measure of entanglement is given by nothing but the 
coherence factor. 

We can explain the different results for one bath and two baths decoherence processes in 
terms of different characters of the Bell states. When the spins interact with separate baths, 
relative orientation of spins is not important because the only difference between the up and 



down configurations is complex conjugation of the coherence factor and it is the modulus 
of the coherence factor which enters the concurrence expression. On the other hand, in the 
single bath case there is no simple relation between the opposite spin terms. In |ei) and |es) 
states, spins are always parallel while in \e 2 ) and \e±) states, they are always antiparallel. 
That is why two groups have different decoherence behaviors. 

Finally, we are going to examine how the Bell inequality is violated in the quantum 
regime and how it is satisfied in the classical domain 15j |. The Bell inequality, or in fact 
inequalities are satisfied if there exists a local realistic theory [la ]. There are a large number 
of Bell inequalities, all resulting from local realistic assumptions, but following Ref. we 
will focus our attention on the quantity 

s = E{e lt e 2 ) - E(e 1 , e' 2 ) + e{9[, 9' 2 ) + e(9[, e 2 ), (20) 

where the correlation function E(8i,8 2 ) is given by, 

E{9 X , 9 2 ) = Tr{c x {ei) ® c 2 (6 2 ) Pc }, (21) 

with 

Ci(0i) = c iz cos 9i + c ix sin 9 { . (22) 
The Bell inequality is violated if \S\ > 2. In calculating whether the inequality is violated, 



es 6i and 9 2 is crucial. It is known that not all entangled states violate 
l^ . That is why #jS must be chosen carefully. In our case, for the Bell 



the choice of the am 
a Bell inequality [18 

states {|ej)}, we can find the corresponding {Si} expressions. To simplify the equations, for 
a given set of angles 9\, 6 2 , 9[, and 9' 2 we will introduce the notation 

A = (cos 9\ cos 9 2 — cos 9\ cos 9' 2 
+ cos 9[ cos 9' 2 + cos 9[ cos 9 2 ) 



(23) 



B = (sin 9\ sin 9 2 — sin 9\ sin 9' 2 
+ sin 9[ sin 9' 2 + sin 9[ sin 9 2 ) 



so that, for two separate baths, 



St = A + B 3?{nr 2 } 
S 2 = -A + B $t{nr*} 
S 3 = A-B ^{nr 2 } 
S 4 = -A- B ^{nr*} 



(24) 
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where r\ and r 2 are again given by Eqn. (|7Jl. and 3?{z} denotes the real part of the complex 
number z. For a single bath very similar expressions hold. In this case 9ft{rir2} and Kjrir^} 
are replaced by 3?{rj^} and 3?{rf 2 }, respectively. 

The angles #1, 6 2 , 6[, and # 2 can take arbitrary values. We are going to pick up a 
particular set for which {Si} are easy to calculate. We will assume that 81 = 0, 9 2 = 7r/4, 
6[ — 7r/2, and # 2 = 37r/4. For this choice of angles, A = B = At t = where all 
decoherence factors are unity, for both two bath and single bath cases, S 2 and 5*3 vanish, 
and therefore they satisfy the Bell inequality |S| < 2. When the system becomes completely 
incoherent so that all coherence factors vanish, again for both two bath and single bath 
cases, we obtain |S 2 | = IS3I = Although there is an increase in |S| values, the inequality 
is still satisfied. In \e\) and |e 4 ) states however, the Bell inequality is violated at t = 0, 
since |Si| = |S 4 | = 2y/2. As the decoherence factors vanish, they both decay to a/2. For 
two baths, I Si I and |S 4 | exhibit the same decay rate. In single bath case, the corresponding 
factors coming from decoherence for |ei) and |e 4 ) states are given by ^{r^} and ^{r^}, 
respectively. As we have discussed above, the two factors decay at different rates. 

In conclusion, using the concurrence and the Bell inequality, we demonstrated that a 
pair of entangled spins show different decoherence behaviors when the spins interact with a 
common spin bath or separate baths. Some entangled states can be more vulnerable than 
others. For example, two entangled electrons in the same quantum dot will have a different 
coherence characteristics than two in separate dots. Recent proposal by Beenakker et al. 
for the creation of entan gled electron-hole pairs might be an interesting system to look for 
such decoherence effects |20j . 
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of the Young Scientist Award Program (MZG/ TUBA-GEBIP/2001-2-9). 
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